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ABSTRACT 



Aims. We present the results of a new spectroscopic and photometric survey of the hot, binary X-ray cluster A115 at z = 0.193, 
containing a radio relic. 

^ , Methods. Our analysis is based on new spectroscopic data obtained at the Telescopic Nazionale Galileo for 115 galaxies and on new 

5_( . photometric data obtained at the Isaac Newton Telescope in a large field. We combine galaxy velocity and position information to 

' select 85 galaxies recognized as cluster members, determine global dynamical properties and detect substructures. 

' Results. We find that Al 15 appears as a well isolated peak in the redshift space, with a global line-of-sight (LOS) velocity dispersion 

CTy = 1362;^|pg km s"'. Our analysis confirms the presence of two structures of cluster-type well recognizable in the plane of the sky 
and shows that they differ of ~ 2000 km s"' in the LOS velocity. The northern, high velocity subcluster (A115N) is likely centred 
on the second brightest cluster galaxy (BCM-A, coincident with radio source 3C28) and the northern X-ray peak. The southern, low 
^ I velocity subcluster (Al 15S) is likely centred on the first brightest cluster galaxy (BCM-B) and the southern X-ray peak. We estimate 

n] , that Al 15S is slightly dynamically more important than Al 15N having CTy = 900- 1 100 km s"' vs. cTv = 750-850 km s"'. Moreover, 

\j ■ we find evidence for two small groups at low velocities. We estimate a global cluster virial mass of 2.2-3.5 x 10'^ h:^^ Mq. 

\l ' Conclusions. Our results agree with a pre-merging scenario where A115N and A115S are colliding with a LOS impact velocity 

Avrf ~ 1600 km s"'. The most likely solution to the two-body problem suggests that the merging axis lies at ~ 20 degrees from the 
1^-^ , plane of the sky and that the cores will cross after ~ 0. 1 Gyr. The radio relic with its largest dimension perpendicular to the merging 

■ axis is likely connected to this merger. 
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1 . Introduction ripheral regions of the cluster Halos and relics are rare sources 

: that appear to be associated with very rich clusters that have un- 

.C^ The evolution of clusters of galaxies as seen in numerical sim- ^^^^^^ ^ ^^^^ ^^^^ suggested by vari- 

ulations IS characterized by the asymmetric accretion ofmass ^^^^^^^^ ^j^^^ ^l^^j^^ halos/relics are related to recent merger 

clumps from surrounding filaments (e.g. Diaferio et al. activity (e.g., Tribble[T993l Burns et al.i994, Fei-etti,1999j. 
Nearby clusters are characterized by a variety of morphologies, 

indicative of different dynamical properties, elongated distribu- The synchrotron radio emission of halos and relics demon- 
tion, traced by several galaxy clumps (e.g., Barrena et al. 120021 strates the existence of large scale cluster magnetic fields, of 
Boschin et al. 120061 Sauvageot et al. 120051 Yuan et al. I2005I ). the order of 0.1-1 fiG, and of widespread relativistic particles 
The presence of substructure is indicative of a cluster in an early of energy density 10 '"^ - 10"'^ erg cm"^. The difficulty in ex- 
phase of the process of dynamical relaxation or of secondary plaining halos/relics arises from the combination of their large 
infall of clumps into already virialized clusters (see Girardi & size and the short synchrotron lifetime of relativistic electrons. 
Biviano |2002| for a review). The expected diffusion velocity of the electron population is on 
An interesting aspect of these investigations is the possible the order of the Alfven speed (~ 100 km s"') making it diffi- 
connection of cluster mergers with the presence of extended, cult for the electrons to diffuse over a megaparsec-scale region 
diff'use radio sources. These sources are large (up to ~ 1 within their radiative lifetime. Therefore, one needs a mecha- 
Mpc), amorphous cluster source s of un certain origin and gener- nism by which the relativistic electron population can be trans- 
ally steep radio spectra (Hanisch [19821 see also and Giovannini ported over large distances in a short time, or a mechanism by 
& Feretti 2002 for a review). They are classified as halos, if they which the local electron population is reaccelerated and the lo- 
are located in the cluster centi-e, or relics, if they appear in the pe- cal magnetic fields are amplified over an extended region. The 

cluster-cluster merger can potentially supply both mechanisms 

Send offprint requests to: R. Barrena, e-mail: rbarrenaOiac . es (e.g., Giovannini et al. ,1993i Burns et al. ,1994i Rottgering et 



2 



Barrena et al.: The dynamical status of Al 15 



al. [19941 see also Feretti et al. |2002l Sarazin |2002l for reviews). 
However, the question is still debated since the diffuse radio 
sources are quite uncommon and only recently we can study 
these phenomena on the basis of a sufficient statistics (a few 
dozens of clusters up to z ~ 0.3, e.g., Giovannini et al. 119991 
see also Giovannini & Feretti |20021 Feretti l2005b . 

Growing evidence of the connection between diffuse ra- 
dio emission and cluster merging is based on X-ray data (e.g., 
Bohringer & Schuecker 120021 Buote l2002b . Studies based on 
a large number of clusters have found a significant relation be- 
tween the radio and the X-ray surface brightness (Govoni et al. 
1200 lal |2001b ) and between the presence of radio-halos/relics 
and irregular and bimodal X-ray surface brightness distribu- 
tion (Schuecker et al. 120011 ). New unprecedented insights into 
merging processes in radio clusters are offered by Chandra 
and XMM-Newton observations (e.g., Markevitch & Vikhlinin 
12001 ; Fuiita etal. 2004). 

Optical data are a powerful way to investigate the presence 
and the dynamics of cluster mergers (e.g., Girardi & Biviano 
120021 ). too. The spatial and kinematical analysis of member 
galaxies allow us to detect and measure the amount of substruc- 
ture, to identify and analyse possible pre-merging clumps or 
merger remnants. This optical information is really complemen- 
tary to X-ray information since galaxies and ICM react on dif- 
ferent time scales during a merger (see numerical simulations by 
Roettiger et al. 1997). Unfortunately, to date optical data is lack- 
ing or poorly exploited and sparse literature concerns some few 
individual clusters. In this context we have caiTied on an inten- 
sive program to study the clusters containing extended, diffuse 
radio emission (Boschin et al. 12004; Boschin et al. 120061 Girardi 
et al. 2006; Barrena et al. 120071 ). In particular, we have conducted 
a study of Abell 115 (hereafter A 11 5) having an extended arc- 
shape reUc (Feretti et al. |1984t Giovannini et al. il987i; Govoni et 
al. l2001bD . 

A115 is a rich cluster (Abell richness class =3; Abell et al. 
11989 ) known in the literature to be characterized by a double 
X-ray peak (A115N and A115S; Forman et al. [TMT] ) and by 
the presence of a strong cooling flow in its northern compo- 
nent (White et al. 119971 ). White et al. ( 119971 ) give the foUow- 
ing values for the X-ray luminosity and temperature: Lx — 14.7 
xlO'*'*, 9.1 xlO"*"* erg s"' (in their cosmology) and Tx - 6.6, 
5.7 keV for the northern and southern components, respectively. 
Slightly smaller temperatures are also found in the literature: 
Tx - 4.9^p g and 5.2 keV for the northern component, and Tx - 
5.2^} J5 and 4.8 keV for the southern component (from Shibata et 
al. |T999i and Gutierrez & Krawczynski 2005 respectively). Both 
ASCA and Chandra data show that the temperature of the ICM is 
highly nonuniform across the cluster suggesting that the merger 
is well underway, but not disturbing the cool subcluster cores 
(Shibata et al. [W99l Gutierrez & Krawczynski 120051 ). From the 
optical point of view. Beers et al. ( 1983) mapped the galaxy dis- 
tribution finding three major clumps of galaxies (A, B, and C in 
Fig. 4a of their paper). Northern and southern clumps A and B 
correspond to the peaks of the X-ray surface brightness distri- 
bution. However, no X-ray emission has been found associated 
with the third eastern clump. 

The northern subcluster contains the very strong radio galaxy 
3C 28 (R iley & Pooley 119751 Feretti et al. [T984l Giovannini et 
al. 119871 ). Its host galaxy is the brightest galaxy of the northern 
subcluster and is classified as "elliptical" by Schombert 119871 
The strong X-ray and radio emissions of this cluster member 
seem to be explained by the presence of a large amount of hot 
gas which is likely to be accreting onto the galaxy. This scenario 
agrees with the observation in the spectrum of this galaxy of the 



Ha-line and other emission lines, likely tracing a cooling flow 
system (Crawford et al. |1999l ). 

The diffuse radio source belongs to A115N, and extends 
from this sub-cluster to the periphery. According to its non- 
central location and its elongated structure, it is classified as a 
cluster relic. However, elongated relics are generally at the clus- 
ter periphery, and with the largest dimension roughly perpen- 
dicular to the cluster radial orientation. So this source is quite 
unusual although several considerations confirm that it should 
be a cluster relic (Govoni et al. l2001bl ). 

To date few spectroscopic data are reported in the literature. 
Beers et al. ( 1983 ) measured the redshift for 19 cluster members 
(see also Zabludoff et al. 119901 for a later re-reduction with a 
new template). In December 2003 we obtained spectra of 115 
galaxies in the cluster region, with the purpose of constraining 
its dynamical status. 

The plan of this paper is the following. In Sect.|2]we describe 
our new photometric and spectroscopic observations. In Sect. [3] 
we provide our results about cluster membership, global proper- 
ties, and substructure. In Sect. |4] we draw our conclusion about 
cluster structure and estimate cluster mass. In Sect. |5] we discuss 
our results and posit a plausible scenario for the dynamical status 
of Al 15. We summarize our results in Sect.|6] 

Throughout this paper, we use - 70 h-iQ km s"' Mpc"' in 
a flat cosmology with n,n - 0.3 and Qa = 0.7. In the adopted 
cosmology, 1' corresponds to ~ 192 kpc at the cluster red- 
shift. 



2. Data 

2.1. Spectroscopy 

Multi-object spectroscopic observations of A115 were caiTied 
out at the TNG telescope in December 2003 during the pro- 
gram of proposal AOT8/CAT-G6. We used DOLORES/MOS 
with the LR-B Grism 1, yielding a dispersion of 187 A/mm, 
and the Loral CCD of 2048 x 2048 pixels (pixel size of 15 
jum). This combination of grating and detector results in disper- 
sions of 2.8 A/pix. We have taken five MOS masks for a total 
of 152 slits. We acquired three exposures of 1800 s for each 
masks. Wavelength calibration was performed using Helium- 
Argon lamps. Reduction of spectroscopic data was carried out 
with IRAFQ package. 

Radial velocities were determined using the cross- 
correlation technique (Tonry & Davis 1979 ) implemented in the 
RVSAO package (developed at the Smithsonian Astrophysical 
Observatory Telescope Data Center). Each spectrum was corre- 
lated against six templates for a variety of galaxy spectral types: 
E, SO, Sa, Sb, Sc, Ir (Kennicutt 7952^. The template produc- 
ing the highest value of % i.e., the parameter given by RVSAO 
and related to the signal-to-noise of the correlation peak, was 
chosen. Moreover, all spectra and their best correlation func- 
tions were examined visually to verify the redshift determina- 
tion. In some cases (IDs 85, 112 and 115, see Table[T]) we took 
the EMSAO redshift as a reliable estimate of the redshift. 

For eleven galaxies we obtained two redshift determinations, 
which are of similar quality. This allow us to obtain a more 
rigorous estimate for the redshift eiTors since the nominal er- 
rors as given by the cross-correlation are known to be smaller 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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than the true errors (e.g., Malumuth et al. 119921 Bardelli et al. 
11994 ', ElUngson & Yee'1994i We fit the first determination vs. 
the second one by using a straight line and considering errors in 
both coordinates (e.g., Press et al. 119921 ). The fitted line agrees 
with the one to one relation, but, when using the nominal cross- 
correlation errors, the small value of ;if^-probability indicates a 
poor fit, suggesting the errors are underestimated. Only when 
nominal errors are multiplied by a ~ 1.3 factor the observed 
scatter can be explained. Therefore, hereafter we assume that 
true errors are larger than nominal cross-correlation errors by a 
factor 1 .3. For the eleven galaxies we used the average of the two 
redshift determinations and the corresponding error 

Our spectroscopic survey consists of 1 15 galaxies taken in a 
field of ~ 15' x20' . 

We also determined the equivalent widths (EW hereafter) of 
the absorption line H6 and the emission line [Oil], in order to 
classify post-starburst and starburst galaxies. We estimated the 
minimum measurable EW of each spectrum as the width of a 
line spanning 2.8 A (our dispersion) in wavelength, with an in- 
tensity three times the rms noise in the adjacent continuum. This 
estimation yields upper measurable limits of ~ 4.5 A in EW. 

We use a conservative approach leading to a sparse spectral 
classification (~ 50% of the sample, see Table [T]i. We follow 
the classification by Dressier et al. (119991 see also Poggianti et 
al.il999). We define "e"-type galaxies those showing active star 
formation as indicated by the presence of an [Oil] and, in partic- 
ular, "e(b)" galaxies when the equivalent width of EW([OII]) 
is < -40 A (likely starburst galaxies); "e(a)"-type galaxies 
those having EW(H<5)>4 A ; "e(c)"-type galaxies those having 
moderate emission lines and EW(H5)<4 A (likely spiral galax- 
ies). We define "k-na" and "a+k"-type galaxies those having 
4 <EW(H5)<8 A and EW(H5)>8 A, respectively, and no emis- 
sion lines (the so called "post-starbust"). Moreover, out of galax- 
ies having the cross-correlation coefficient 'R>1- correspond- 
ing to S /N > 10 as obtained using A773 data (Barrena et al. 

I2007I I - we define "k"-type galaxies those having EW(Hi5)<3 A 
and no emission lines (likely passive galaxies). We classify 62 
galaxies finding 38 "passive" galaxies, 24 "active" galaxies (i.e., 
14 "k-i-a"/"a-i-k" and 10 "e" galaxies). 

2.2. Photometry 

As far as photometry is concerned, our observations were car- 
ried out with the Wide Field Camera (WFC), mounted at the 
prime focus of the 2.5m INT telescope (located at Roque de los 
Muchachos observatory. La Palma, Spain). We observed A115 
in 1 8 December 2004 in photometric conditions with a seeing of 
about 2". 

The WFC consists of a 4 chip mosaic covering a 
30'x30' field of view, with only a 20% marginally vignetted 
area. We took 10 exposures of 720 s in Bu and 360 s in Ru 
Harris filters (a total of 7200 s and 3600 s in each band) de- 
veloping a dithering pattern of ten positions. This observing 
mode allowed us to build a master "supersky" image that was 
used to correct our images for vignetting and fringing patterns 
(Gullixson |1992| l. In addition, the dithering helped us to clean 
cosmic rays and avoid gaps between CCD chips. The complete 
reduction process (including flat fielding, bias subtraction and 
bad-columns elimination) yielded a final co-added image where 
the variation of the sky was lower than 0.4% in the whole frame. 
Another effect associated with the wide field frames is the dis- 
tortion of the field. In order to match the photometry of several 
filters (in our case, only Z?h and Ru), a good astrometric solu- 



tion taking into account these distortions is needed. Using IRAF 
tasks and taking as reference the USNO Bl.O catalog we were 
able to find an accurate astrometric solution (rms~ 0.5") across 
the full frame. The photometric calibration was performed using 
Landolt standard fields achieved during the observation. 

We finally identified galaxies in our Bu and Ru images and 
measured their magnitudes with the SExtractor package (Bertin 
& Arnouts 1996) and AUTOMAG procedure. In few cases, (e.g., 
close companion galaxies, galaxies close to defects of CCD), 
the standard SExtractor photometric procedure failed. In these 
cases we computed magnitudes by hand. This method consists 
in assuming a galaxy profile of a typical elliptical and scale it 
to the maximum observed value. The integration of this profile 
give us an estimate of the magnitude. The idea of this method 
is similar to the PSF photometry, but assuming a galaxy profile, 
more appropriate in this case. 

We transformed all magnitudes into the Johnson- 
Cousins system (Johnson & Morgan 119531 Cousins 
[T976] l. We used B ^ Bu + 0.13 and R = Ru, 
as derived from the Harris filter characterization 
(http://www.ast.cam.ac.uk/~wfcsur/technical/photom/colours/) 
and assuming a B-V ~ 1.0 for E-type galaxies (Poggianti 
1997 ). As a final step, we estimated and corrected the galactic 
extinction Ab ~ 0.25, Ar ~ 0.15 from Burstein & Heiles (|1982| l 
reddening maps. We estimated that our photometric sample is 
complete down to B = 19.5 (21.0) and R = 22.0 (23.0) for 

5 /N - 5 (3) within the observed field. 

We assigned B and R magnitudes to the whole spectro- 
scopic sample. We measured redshift for galaxies down to R ~ 
20.5 mag, but we are complete to 60% down to R =18 mag 
[within a region of 13' x 20 ' around R.A.=00*'56"02 ?0, 
Dec.=-H26°23' 00" (J2000.0)]. 

Table [T] lists the velocity catalogue (see also Fig. [T] and 
Fig. |2|i: identification number of each galaxy, ID (Col. 1); ID 
code following the lAU nomenclature (Col. 2); right ascen- 
sion and declination, a and 6 (J2000, Col. 3); B and R magni- 
tudes (Cols. 4 and 5, respectively); heliocentric radial velocities, 
V = czQ with errors, Av (Cols. 6 and 7, respectively); spectral 
classification SC (Col. 8). 

3. Analysis & Results 

3. 1 . Member selection and global properties 

To select cluster members out of the 115 galaxies having red- 
shifts, we follow the two steps procedure. First, we perform the 
adaptive-kernel method (hereafter DEDICA, Pisani 119931 and 
1996; see also Fadda et al. 1996; Girardi et al. [T9961 Girardi 

6 Mezzetti .2001 ). We find the significant peaks in the velocity 
distribution >99% c.l.. This procedure detects A115 as a one- 
peak structure at z ~ 0. 1937 populated by 88 galaxies considered 
as candidate cluster members (see Fig.O. Out of non-member 
galaxies, 12 and 15 are foreground and background galaxies, re- 
spectively. 

The projected clustercentric distance vs. the rest-frame ve- 
locity of the 88 galaxies is shown in Fig. |4] where we also 
show the three brightest galaxies in our sample, each of them 
corresponding to a different galaxy density clumps as found 
by Beers et al. ( IT983b . i.e. IDs 54, 21(=3C28) and 104 (here- 
after BCM-B, BCM-A, and BCM-C) coiTesponding to B 
(A115S), A (A115N), and C clumps, respectively. The posi- 
tion of BCM-B is very close to that of the southern X-ray peak 
[R.A.=00'>55™58 ?81, Dec.=-H26° 19' 58 '.'9 (J2000.0) as recov- 
ered from X-ray Chandra data, see Fig.|2|. BCM-A=3C28 is co- 
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Fig. 3. Redshift galaxy distribution. The solid line histogram 
refers to galaxies assigned to the cluster according to the 
DEDICA reconstruction method. 



incident with the northern X-ray peak and is very close to the po- 
sition of the X-ray centroid when masking the X-ray emission 
of the point source [R.A.^00^^53 ?2, Dec. =+26° 24' 59" 
(J2000.0) by Govoni et al. 1 200 Ibl . We also consider the bright 
galaxy ID 81 (hereafter BCM-D, see Sect. [33]l located in the 
middle of the above three galaxies. 

All the galaxies assigned to the A115 peak are analysed in 
a second step by applying the "shifting gapper" technique by 
Fadda et al. ( 1996 ), which uses the combination of position and 
velocity information. This procedure rejects galaxies that are too 
far in velocity from the main body of galaxies and within a fixed 
bin that shifts along the distance from the cluster centre. The pro- 
cedure is iterated until the number of cluster members converges 
to a stable value. Following Fadda et al. ( 1996) we use a gap of 
1000 km s"' - in the cluster rest-frame - and a bin of 0.6 /i^q 
Mpc, or large enough to include 15 galaxies. In the case of the 
A115 complex, which is not an individual system, we fix alter- 
natively BCM-A and -B as cluster centres. Fig.|4]shows the re- 
sults. The very high velocity galaxy (ID 91) is a clear interloper 
far more than 3000 km s"' from other galaxies. The other two 
high velocity galaxies (IDs 96 and 100), that are close enough in 
2D and far from the high velocity BCM-A galaxy (see Fig. [TJ, 
are rejected in both our analyses. The conclusion about the low 
velocity tail is less clear also because few of these galaxies lie 
in the eastern region of the cluster and might be associated to 
the C clump (or to the D clump, see Sect. 13.31 ). We decide to de- 
fine 85 likely cluster members rejecting the three highest veloc- 
ity galaxies (basic sample). In some analyses we also consider 
an alternative sample of 80 galaxies - Sample2 - also rejecting 
the low velocity galaxies indicated in Fig. H] (the crosses in the 
middle panel). 

By applying the bi-weight estimator to cluster members 
(Beers et al. 1990), we compute a mean cluster redshift of 
<z) = 0.1929+ 0.0005, i.e. (v) = 57817+148 km s"'. We es- 
timate the LOS velocity dispersion, ctv, by using the bi-weight 
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Fig. 4. Projected clustercentric distance vs. rest-frame velocity 
Vrf =(v-<v>)/(l + z) for the 88 galaxies in the main peak (see 
Fig. [3]) showing galaxies detected as interlopers by our "shift- 
ing gapper" procedure (stars and crosses) when using BCM-A 
and -B as alternative cluster centres (top and bottom panels, re- 
spectively). We reject the three galaxies with highest velocities 
(stars) to build our basic sample of 85 galaxies and all galax- 
ies indicated by crosses and stars in the top panel to build our 
Sample2. The three brightest cluster members each one corre- 
sponding to Beers et al. ( 119831 ) clumps (BCM-A, -B and -C) 
and a fourth bright galaxy (BCM-D) are indicated by squares. 



estimator and applying the cosmological coiTection and the stan- 
dard coiTection for velocity eiTors (Danese et al. |1980l ). We ob- 
tain (Tv = 1362^}^^ km s"', where eiTors are estimated through 
a bootstrap technique. Consistent results are found for Sample2 
(see Table|2|). 

Hereafter, for practical reasons, we consider as centre of the 
whole A115 complex the position of the bi-weight centre ob- 
tained using bi-weight mean estimators for R.A. and Dec. sepa- 
rately [R.A.=00'>56™01 ?31, Dec.=+26°22'26 :'6 (J2000.0)]. 



3.2. Velocity structure 

We analyse the velocity distribution to look for possible devia- 
tions from Gaussianity that could provide important signatures 
of complex dynamics. For the following tests the null hypothesis 
is that the velocity distribution is a single Gaussian. 

We estimate three shape estimators, i.e. the kurtosis, the 
skewness, and the scaled tail index (see, e.g.. Beers et al. |1991| ). 
The value of the skewness (-0.462) shows evidence that the 
velocity distribution differs from a Gaussian at the 95 - 99% 
c.l. (see Table 2 of Bir d & Beers 119931 ). Moreover, the W-test 
(Shapiro & Wilk |1965] ) marginally rejects the null hypothesis of 
a Gaussian parent distribution at the 92% c.l.. 

Then we investigate the presence of gaps in the distribution. 
A weighted gap in the space of the ordered velocities is defined 
as the difference between two contiguous velocities, weighted 
by the location of these velocities with respect to the middle of 
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Table 2. Results of the kinematical analysis. 
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We use the biweigth and the gapper estimators by Beers et al. 
(1990) for samples with Ng > 15 and with Ng < 15 galaxies, re- 
spectively (see also Girardi et al. ll993l l. 

Notice that BCM-B resembles characteristics of a passive galaxy, 
but it does not appears in our classification since its spectra has H 
slightly below the threshold value of 7. 
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Fig. 5. Top panel: velocity distribution of the 85 fiducial clus- 
ter members. Arrows correspond to the bright galaxies BCM-A, 
-B, -C and -D). Bottom panel: stripe density plot where the ar- 
rows indicate the positions of the significant gaps. The gap at the 
lower velocity has a normalized size =2.41, the other =2.25). 



the data. We obtain values for these gaps relative to their av- 
erage size, precisely the midmean of the weighted-gap distri- 
bution. We look for normalized gaps larger than 2.25 since in 
random draws of a Gaussian distribution they arise at most in 
about 3% of the cases, independent of the sample size (Wainer 
and Schacht 1978; see also Beers et al. 1991). Two significant 
gaps in the ordered velocity dataset are detected individuating 
two groups in the low velocity tail of the velocity distribution, of 
eight and six galaxies (see WGAPl and WGAP2 in Table|2]and 
Fig.E). 

We use the results of the gap analysis to determine the first 
guess when using the Kaye's mixture model (KMM) to find a 
possible group partition of the velocity distribution (as imple- 
mented by Ashman et al. 1994). The KMM algorithm fits an 
user-specified number of Gaussian distributions to a dataset and 
assesses the improvement of that fit over a single Gaussian. In 
addition, it provides the maximum-likelihood estimate of the 
unknown n-mode Gaussians and an assignment of objects into 
groups. We do not find any two- or three-groups partition which 
is a significant better descriptor of the velocity distribution with 
respect to a single Gaussian. 



3.3. 2D galaxy distribution 

When applying the DEDICA method to the 2D distribution of 
A115 galaxy members we find three significant peaks. The po- 
sition of the highest peak is close to the location of the BCM- 
B and of the southern peak of X-ray emission (XS in Fig|6]l. 
Another peak is close to the location of the BCM-A, of the 
northern peak of X-ray emission (XN) and of the X-ray cen- 
troid when masking the X-ray emission of the point source (X, 
Govoni et al. 2001b). When dividing the sample in bright and 
faint galaxies - using the median magnitude value /?=18.32 - we 
find that the 2D distributions of the two samples are different at 



the 98.5% c.l. according to the two-dimensional Kolmogorov- 
Smirnov test (hereafter 2DKS-test; see Fasano & Franceschini 
119871 as implemented by Press et al. il992i) . In fact, the faint 
galaxies sample shows both peaks A and B, while the bright 
galaxies sample only shows the peak B. 

Our spectroscopic data do not cover the entire cluster field 
and suffer from magnitude incompleteness. To overcome these 
limits we recover our photometric catalogue selecting likely 
members on the basis of the colour-magnitude relation (here- 
after CMR), which indicates the early-type galaxy locus. To de- 
termine CMR we fix the slope according to Lopez-Cruz et al. 
(120041 see their Fig. 3) and apply the two-sigma-clipping fit- 
ting procedure to the cluster members obtaining B-R = 3.61 1 - 
0.069 xR (see Fig.|7]i. Out of our photometric catalogue we con- 
sider galaxies (objects with SExtractor stellar index < 0.9) ly- 
ing within 0.25 mag of the CMR. To avoid contamination by 
field galaxies we do not show results for galaxies fainter than 21 
mag (in /?-band). The contour map for 369 likely cluster mem- 
bers having R < 21 shows again the two peaks in coiTespon- 
dence of BCM-A and BCM-B and also a peak in coiTespon- 
dence of BCM-C. Finally, a less dense peak lies in the middle 
of the above three peaks and corresponds to the luminous galaxy 
BCM-D at a very low velocity (see Fig. |8]l. Similar results are 
found analysing the 268 galaxies with R < 20. The analysis of 
136 galaxies with R < 19 shows as very significant only the 
southern peak close to BCM-B in agreement with the results 
coming from the spectroscopic sample (see above). 



3.4. Position-velocity correlations 

The existence of correlations between positions and velocities of 
cluster galaxies is a footprint of real substructures. Here we use 
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Fig. 6. Spatial distribution on the sky of spectroscopically con- 
firmed cluster members and the relative isodensity contour map. 
The three brightest galaxies BCM-A -B, and -C corresponding 
to Beers et al. Q983 ) groups, a fourth bright galaxy (BCM-D) 
and the X-ray peaks are indicated by large squares, too. The plot 
is centred on the cluster centre defined as the bi-weight centre 
(see text). 
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Fig. 7. B-R vs. R diagram for galaxies with available spec- 
troscopy is shown by small dots and crosses (cluster and field 
members, respectively). Large solid dots indicate luminous 
galaxies BCM-A, -B, -C and -D. The solid line gives the best- 
fit colour-magnitude relation as determined on member galax- 
ies; the dashed lines are drawn at +0.25 mag from the CMR. 
Large open symbols indicates classified galaxies: circles ("k"), 
squares ("k-na"), rombs ("a-nk"), triangles ("e(c)"), and stars 
("e(a)"/"e(b)"). 
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Fig. 8. Spatial distribution on the sky and relative isodensity con- 
tour map of 369 likely cluster members (according to the colour- 
magnitude relation) with R < 21, obtained with the DEDICA 
method. The three brightest galaxies BCM-A -B, and -C are 
indicated by large circles and the X-ray peaks are indicated 
by large squares. It is also shown the position of the luminous 
galaxy BCM-D. The plot is centred on the cluster centre. 



different approaches to analyse the structure of A115 combining 
velocity and position information. 

The cluster velocity field may be influenced by the presence 
of internal substructures. To investigate the velocity field of the 
All 5 complex we divide galaxies in a low and a high velocity 
samples by using the median cluster velocity and check the dif- 
ference between the two distributions of galaxy positions. Fig.|9] 
shows that low and high velocity galaxies are segregated roughly 
in the E-W direction. The two distributions are different at the 
99.5% c.l. according to the 2DKS-test. In order to estimate the 
direction of the velocity gradient we perform a multiple linear re- 
gression fit to the observed velocities with respect to the galaxy 
positions in the plane of the sky (see also den Hartog & Katgert 
1996 . Girardi et al. 1 19961 1. We find a position angle on the celes- 
tial sphere of PA = 269!j^ de grees (measured counter-clock- 



wise from north), i.e. higher velocity galaxies lie in the western 
region of the cluster (see Fig. |9l). To assess the significance of 
this velocity gradient we perform 1000 Monte Carlo simulations 
by randomly shuffling the galaxy velocities and for each sim- 
ulation we determine the coefficient of multiple determination 
(RC^, see e.g., NAG Fortran Workstation Handbook 1986^. We 
define the significance of the velocity gradient as the fraction of 
times in which the RC^ of the simulated data is smaller than the 
observed RC^. We find that the velocity gradient is marginally 
significant at the 90% c.l.. Similar results are obtained for the 
Sample2 {PA = 267+}^ degrees at the 91% c.l.). 

We combine galaxy velocity and position information to 
compute the A-statistics devised by Dressier & Schectman 
(1988J . This test is sensitive to spatially compact subsystems 
that have either an average velocity that differs from the cluster 
mean, or a velocity dispersion that differs from the global one, 
or both. We find A = 121 for the value of the parameter which 
gives the cumulative deviation. This value is an indication of 
substructure, significant at the 96% c.l., as assessed computing 
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Fig. 9. Spatial distribution on the sky of 85 cluster members. 
Open and solid circles indicate low and high velocity galaxies, 
the cross the galaxy with median velocity. The larger the sym- 
bol, the smaller is the radial velocity. The solid and faint lines 
indicate the position angle of the cluster gradient and relative er- 
rors, respectively. The faint big squares indicate the 14 galaxies 
belonging to WGAPl and WGAP2 (see Sect.[l2]l- The plot is 
centred on the cluster centre. 



1000 Monte Carlo simulations, randomly shuffling the galaxy 
velocities. Fig. [TOl shows the distribution on the sky of all galax- 
ies, each marked by a circle: the larger the circle, the larger the 
deviation i5, of the local parameters from the global cluster pa- 
rameters, i.e. the higher the evidence for substructure. This fig- 
ure provides information on the positions of substructures: one 
in the eastern region corresponding to the clump C and one in 
the northern region corresponding to the clump A. Similar re- 
sults obtained for Sample2 (A = 1 10 and a c.l. of 95%), but the 
eastern substructure is no longer so obvious in the relative plot. 

To obtain further information, we resort to the technique de- 
veloped by Biviano et al. (2002^, who used the individual <5,- 
values of the Dressier & Schectman method. The critical point is 
to determine the value of <5/ that optimally indicates galaxies be- 
longing to substructure. To this aim we consider the i5,— values of 
all 1000 Monte Carlo simulations already used to determine the 
significance of the substructure (see above). The resulting distri- 
bution of 6i is compared to the observed one finding a difference 
of P > 99.99% c.l. according to the KS-test. The "simulated" 
distribution is normalized to produce the observed number of 
galaxies and compared to the observed distribution in Fig. [TT] 
the latter shows a tail at large values. The tail with 5 > 2 is pop- 
ulated by galaxies that presumably are in substructures. For the 
selection of galaxies within substructures we choose the value of 
Slim - 2.2, since the galaxies with <5/ > dum are well separated in 
the sky (see Fig.fTOll and assign three and eleven galaxies to the 
clumps C and A, respectively (hereafter DS-C and DS-A). The 
northern structure is populated by high velocity galaxies, while 
the poor statistics prevent us to obtain firm conclusions about the 
eastern structure (see Table |2]i. 



Fig. 10. Spatial distribution of cluster members, each marked by 
a circle: the larger the circle, the larger is the deviation (5, of the 
local parameters from the global cluster parameters, i.e. there 
is more evidence for substructure (according to the Dressier & 
Schectman test, see text). The boldface circles indicate those 
with di > 2.2 (see text). The plot is centred on the cluster centre. 



The Gaussian model for the 2D galaxy distribution is poorly 
supported by theoretical and/or empirical arguments and, how- 
ever, our galaxy catalogue is not complete down to a magnitude 
limit. However, since the 3D diagnostics is in general the most 
sensitive indicators of the presence of substructure (e.g., Pinkney 
et al. 1996 ), we apply the 3D version of the KMM software using 
simultaneously galaxy velocity and positions. We use the galaxy 
assignment given by Dressler-Schectman method to determine 
the first guess when fitting three groups. The algorithm fits a 
three-groups partition at the 97% c.l. according to the likelihood 
ratio test (hereafter KMMl, KMM2, KMM3 groups from low to 
high mean velocities). The results for the three groups are shown 
in Table |2] and Fig. [12] KMMl group is sparse in the sky, but 
well distinct in velocity from the other two groups. Several of 
its galaxies were already detected by the weighted gap analysis 
as belonging to WGAPl and WGAP2 (see Sect.[3l2]i. Moreover, 
KMMl contains BCM-D. KMM2 and KMM3 groups ai-e well 
distinct in the sky. KMM2 contains both BCM-B and BCM-C, 
while KMM3 contains BCM-A. 

3.5. Kinematics of A115N and A115S 

The spatial agreement between the two brightest cluster mem- 
bers (BCM-A and -B) and the peaks of X-ray emission as well 
as the high density of galaxies around BCM-A and -B prompts 
us to analyse the profiles of the mean velocity and the veloc- 
ity dispersion of galaxy systems surrounding these two galaxies 
(see Figs. [13] and [141 respectively). This allows an independent 
analysis of the individual galaxy clumps. Although an increase 
in the velocity-dispersion profile in the cluster central regions 
might be due to dynamical friction and galaxy merging (e.g., 
Menci & Fusco-Femiano 1 1 9961 Girardi et al. 119981 Biviano & 
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Fig. 11. The distribution of <5, deviations of the Dressler- 
Schectman analysis for the 85 member galaxies. The solid line 
represents the observations, the dashed line the distribution for 
the galaxies of simulated clusters, normalized to the observed 
number. 



Fig. 12. Spatial distribution on the sky of the 85 member galax- 
ies. Solid circles, open circles and triangles indicate KMMl, 
KMM2, and KMM3. The large faint squares indicate the po- 
sition of the brightest cluster members BCM-A, -B, -C and -D. 
The plot is centred on the cluster centre. 



Katgert |20041 l, in the case of A115 it might be simply induced 
by the contamination of the galaxies of a secondary clump (e.g., 
Girardi et al.'1996i Girardi et al. 2005 ). The latter hypothesis can 
be investigated by looking at the behaviour of the mean velocity 
profile. 

Since from the above sections we know the presence of one 
or likely a few low velocity groups, we analyse the Sample2 to 
avoid, at least partially, the possible contamination. We also con- 
sider the results obtained rejecting all galaxies belonging to low 
velocity WGAPl and WGAP2. Figs. [n]and[T3show velocity- 
dispersion and mean-velocity profiles, as well as regions not 
likely to be contaminated by other galaxy systems and thus reli- 
able for kinematical analysis. Detailed results of this analysis are 
included in Table|2]where the "uncontaminated" galaxy systems 
are named as CORE-A and CORE-B. 

Fig. [T3] shows how the integral velocity-dispersion slightly 
increases with the distance from the BCM-A. Simultaneously, 
the mean velocity shows a continuous decline from high values 
~ 59000 km s"' suggesting a strong contamination of galaxies 
from structures connected with BCM-B, -C and -D, all having 
lower mean velocities. In a conservative view we consider the 
likely uncontaminated region within 0.25 /z^q Mpc, where we 
find cTv - 750 km s"'. Fig. [T4l shows an enough robust mean 
velocity and a sharp increase of the integral velocity-dispersion 
with the distance from the BCM-B up to a peak value ~ 1450 
km s"' at ~ 0.3 Mpc. We interpret these features as the con- 
tamination of the structures connected with BCM-A, -C and -D 
with higher and lower mean velocities. Their combination does 
not affect the mean velocity but strongly increases the veloc- 
ity dispersion. In fact, the peak value of the velocity dispersion 
goes down from 1450 to 1200 km s"' when rejecting low ve- 
locity galaxies of WGAPl and WGAP2. We then consider the 



likely uncontaminated region within 0.25 h^^ Mpc, where we 
findo-v ^900kms-i. 



4. 3D structure and virial mass of A115 

On the basis of the above section we conclude that A115 is 
formed by two subclusters well distinct in the sky and centred 
around BCM-A and -B, hereafter Al 15N and Al 15S, with the 
addition of several low velocity galaxies, above all in the eastern 
cluster region, likely organised in small groups (see Fig. |9] and 
the PA of the velocity gradient). 

As for the low velocity groups we detect two 2D galaxy con- 
centrations around the two bright, low velocity galaxies BCM-C 
and -D. Moreover, we find two possible groups in the low veloc- 
ity tail of the velocity distribution (WGAPl containing BCM-D 
and WGAP2). The Dressler-Schectman test detects a substruc- 
ture around the BCM-C, too. The presence of a group around 
BCM-D finds some support in the faint X-ray emission shown 
there (see Gutierrez & Krawcz ynski 120051 and Fig.|2]i. That we 
find evidence of a few groups rather than of an individual mas- 
sive cluster agrees with the absence of a third X-ray luminous 
peak. 

As for Al 15N, we find comparable kinematical properties in 
the analysis of DS-A and CORE-A, i.e. <v) > 59000 km s ' and 
(Tv ~ 800 km s"'. Instead KMM3 is contaminated by the pres- 
ence of low velocity galaxies, e.g. IDs 10 and 19 close to BCM- 
A. As for A115S, we find (v> ~ 57500 km s"' and cr^ ~ 900- 
1000 km s ' (see KMM2 and CORE-B). The velocity disper- 
sion of the two subclusters are roughly comparable to their aver- 
age X-ray temperature as listed in the literature and transformed 
in assuming the density-energy equipartition between gas 
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Fig. 13. Kinematical profiles of the northern subcluster A115N 
obtained assuming BCM-A as centre. The vertical line indicates 
the region likely not contaminated from other galaxy clumps 
(see text). Top panel: rest-frame velocity vs. projected distance 
from the subcluster centre (BCM-A): crosses indicate the galax- 
ies belonging to WGAPl and WGAP2. Differential (big circles) 
and integral (small points) mean velocity and LOS velocity- 
dispersion profiles are shown in middle and bottom panels, re- 
spectively. For the differential profiles we plot the values for 
seven annuli from the centre of the subcluster, each of 0.25 h^^^ 
Mpc. For the integral profiles, the mean and dispersion at a given 
(projected) radius from the subcluster-centre is estimated by 
considering all galaxies within that radius - the first value com- 
puted on the five galaxies closest to the centre. The error bands at 
the 68% c.l. are also shown. In the bottom panel, the horizontal 
line represents the range of X-ray temperatures as given in the 
literature for A115N (see Sect. [U transformed in o-y assuming 
the density-energy equipartition between gas and galaxies, i.e. 
Apec = 1 (see text). 



and galaxies, i.e. ySspec = 10 (see Figs. [T3] and [T4li. The two sub- 
clusters differ by ~ 2000 km s"' in the LOS velocity, i.e. almost 
three times more than what found by Beers et al. (11983t with 
only 19 cluster members. In agreement with Beers et al. (I1983I I, 
we find that A115S is dynamically somewhat more important 
than A115N, while the contrary is found by X-ray data (Jones 
& Forman 119991 White et al. 119971 GutieiTez & Ki-awczynski 
120051 but see Shibata et al. ll999l l. As already suggested by Beers 
et al. (I1983I I the presence of 3C28 might affects the X-ray re- 
sults overestimating the X-ray luminosity and the temperature 
of A115N. 

Although A115 is likely in a phase of interaction (see the 
following section), the two main galaxy subclusters are still well 
detectable and roughly well coincident with the X-ray peaks. 
Thus we assume that each subcluster is in dynamical equilib- 
rium to compute virial quantities. Hereafter we assume cTy - 



^ Apec = c^likT I^Mp) with = 0.58 the mean molecular weight and 
nip the proton mass. 



Fig. 14. The same as in Fig. [13] but referring to the southern sub- 
cluster A115S centred around BCM-B. The dashed line in the 
bottom panel gives the integral velocity dispersion profile when 
rejecting the low velocity galaxies of WGAPl and WGAP2. In 
the bottom panel, the horizontal line represents the range of X- 
ray temperatures as given in the literature for Al 15S (see Sect.lTJ 
transformed in cTv assuming the density-energy equipartition be- 
tween gas and galaxies, i.e. ySspec = 1 (see text). 



750-850 km s~i and 900-1000 km s"' for Al 15N and Al 15S, 
respectively. 

Following the prescriptions of Girardi & Mezzetti ( I2001I I. 
we assume for the radius of the quasi-virialized region Ryir n = 

0. 17 X crv///(z) = 2.2-2.5 /i",; Mpc and Rvi,-,s = 2.6-2.9 h'^ 
Mpc for Al 15N and Al 15S, respectively - see their eq. 1 after 
introducing the scaling with H{z) (see also eq. 8 of Carlberg et 
al. 1 19971 for R2oo)- Therefore, our spectroscopic catalogue sam- 
ples most of the virialized region. 

One can compute the mass using the virial theorem (Limber 
& Mathews fl 9601 see also, e.g., Girardi et al l 19981 1 under the 
assumption that mass follows galaxy distribution: M - Msvir - 
SPT, where Msvh- - 37r/2 • cr^Rpy/G is the standard virial mass, 
Rpv a projected radius (equal two times the harmonic radius), 
and SPT is the surface pressure term correction (The & White 
I1986I I. The value of Rpv depends on the size of the considered re- 
gion, so that the computed mass increases (but not linearly) with 
the increasing considered region. Since the two subclusters are 
interacting and we do not cover the whole virialized region we 
use an alternative estimate which was shown be good when Rpv 
is computed within Ryir (see eq. 13 of Girardi et al. |1998l l. This 
alternative estimate is based on the knowledge of the galaxy dis- 
tribution and, in particular, a galaxy King-like distribution with 
parameters typical of nearby/medium-redshift clusters: a core 
radius Rcore = 1/20 x Ryir and a slope-parameter jS^t - 0.8, 

1. e. the volume galaxy density at large radii as r^^l^^'^ - r^^-'^ 
(Girardi & Mezzetti '2001). For the whole virialized region we 
obtain Rpv,n - 1.6-1.8 /z^q Mpc and Rpv,s =1.9-2.2 Mpc. 
As for the SPT correction, we assume a 20% computed com- 
bining data on many clusters one (e.g., Carlberg et al. 119971 
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Girardi et al. I1998I I. This leads to virial masses Mn(< Rvir.N - 
2.2 - 2.5 h-^Mpc) = 8.0-1 1.6 xlO'"* h'^ Mq and Ms(< Rvir.s = 
2.6 - 2.9 /!7(}Mpc) = 13.8 -18.9 xlO'"* h'^ for the two 
subclusters. 

To compare our results with the estimate recovered from 
a X-ray surface brightness deprojection analysis (White et al. 
11997 ) we assume that each subcluster is described by a King- 
like mass distribution (see above) or, alternatively, a NFW pro- 
file where the mass-dependent concentration parameter is taken 
from Navarro et al. ( 1 19971 1 and rescaled by the factor 1 + z 
(Bullock et al. 120011 Dolag et al. 2004). We obtain Mn(< 
Rout,N = 0.235 /i^o'Mpc) = (4.2 - 9.9) xlO'^ /j^o' M^and 
Ms(< Rout,s = 0.157 h-^Mpc) = (2.5 - 5.5) xlO'^ h'^ Mq. 
The first estimate is somewhat smaller and the second one is 
in agreement with those found by White et al. ( 119971 see their 
Table 3 where Mgi-av,N(< Rout,N) - 14.0 xlO'^ h'^ Mq and 
A^grav,s(< Rout,s) - 5.0 xlO'^ Mq in our cosmology). 

As for the mass of the whole system, the contribution of 
the low velocity groups is of minor importance since they likely 
have low velocity dispersion and the virial mass scales with crl. 
E.g., we estimate that M ~ 2.4 xlO''* /i^,,' Mq for a ctv ~ 500 
km s ' group (see WGAPl in Table |2|. Considering the possi- 
ble presence of, at most, two of these cr^ ~ 500 km s"' groups, 
a reliable mass estimate of the whole system is then M=2.2- 
3.5 xlO'^ li^Q Mq, in agreement with rich clusters reported 
in the literature (e.g., Girardi et al. '1998'; Girardi & Mezzetti 
120011 ). A smaller mass estimate is given by Govoni et al. 1 200 Ibl 
M(< R = 2.3 h^y'Mpc) = 0.6-1.2 xlO'^ /i^d centred on the 
X-ray centroid in Al 15N], but notice that it refers to a smaller 
cluster region, i.e. likely excluding a large part of Al 15S. 

5. Merging scenario 

Since Al 15N and Al 15S are well detectable and optical and X- 
ray data indicate a very similar location we are likely looking at 
the cluster prior to merging. However, A115N and A115S sub- 
clusters are already starting to interact, as suggested by several 
pieces of evidence: the slight displacement between peaks of gas 
distribution and of galaxy distribution (see our Fig. [8]); the pres- 
ence of an hot region likely due to the interaction (Shibata et 
al. 119991 GutieiTez & Krawczvnski 12005]) . Very noticeably, the 
largest dimension of the radio relic is somewhat perpendicular 
to the axis connecting Al 15N and Al 15S in agreement with be- 
ing originated by shock waves connected to the ongoing merger 
(e.g., Ensslin & Bruggen l2002l ). 

When the merging scenario is assumed to explain the pres- 
ence of the hot region located between A115N and A115S, a 
relative colliding velocity is necessary to heat up the ICM tem- 
perature kT from ~ 4 to ~ 9 keV (see Gutierrez & Krawczynski 
12005 ). Assuming that the two subclusters are to cause a head- 
on collision and that their kinetic energies are completely con- 
verted to thermal energy, the necessary value of the colliding 
velocity is Vcou = (3/tAr/jUOTp)'/2 (^gg Shibata et al. |T999l ). We 
find VcoU ~ 1600 km s ' in good agreement with the observed 
relative LOS velocity in the rest frame Vi - 1646 km s"' as re- 
covered from CORE-A and CORE-B in the cluster rest frame, 
i.e. (<v>core-a-<v>core-b)/(1 + z)). 

5.1. Bimodal model 

Here we investigate the relative dynamics of Al 15N and Al 15S 
using different analytic approaches which are based on an energy 
integral formalism in the framework of locally flat spacetime and 



Newtonian gravity (e.g.. Beers et al. 119821 ). The values of the rel- 
evant observable quantities for the two-clumps system are: the 
relative LOS velocity in the rest frame, Vi = 1646 km s"' (as re- 
covered from CORE-A and CORE-B); the projected linear dis- 
tance between the two clumps, D = 0.89 /i^^ Mpc (as recovered 
from BCM-A and -B cluster rest frame); the mass of the sys- 
tem obtained by adding the masses of the two subclusters each 
within its virial radius, logMsys = 15.4107+qq^34 (see Sect.©. 

First, we consider the Newtonian criterion for gravita- 
tional binding stated in terms of the observables as V^D < 
2GMsys sin^ a cos (T, where a is the projection angle between 
the plane of the sky and the line connecting the centres of two 
clumps. The faint curve in Fig. [TSl separates the bound and un- 
bound regions according to the Newtonian criterion (above and 
below the curve, respectively). Considering the value of Msys, 
the A115N-I-S system is bound between 20° and 84°; the cor- 
responding probability, computed considering the solid angles 

(i.e., JjQo cos a da), is 65%. We also consider the implemented 

criterion V^D < 2GM sin^ cos o-d, which introduces different 
angles and for projection of distance and velocity, not as- 
suming strictly radial motion between the clumps (Hughes et al. 
1995J . We obtain a binding probability of 60%. 

Then, we apply the analytical two-body model introduced 
by Beers et al. (119821 ) and Thompson (119821 see also Lubin et 
al. il998i for a recent application). This model assumes radial or- 
bits for the clumps with no shear or net rotation of the system. 
Furthermore, the clumps are assumed to start their evolution at 
time to - with separation = 0, and are moving apart or com- 
ing together for the first time in their history; i.e. we are assum- 
ing that we are seeing the cluster prior to merging (at the time 
t=11.106 Gyr at the cluster redshift). The bimodal model solu- 
tion gives the total system mass Msys as a function of a (e.g., 
Gregory & Thompson 1984) . Fig. [15] compares the bimodal- 
model solutions with the observed mass of the system, which is 
the most uncertain observational parameter The present bound 
outgoing solutions (i.e. expanding), BO, are clearly inconsistent 
with the observed mass. The possible solutions span these cases: 
the bound and present incoming solution (i.e. collapsing), Bla 
and Bib, and the unbound-outgoing solution, UO. For the in- 
coming case there are two solutions because of the ambiguity 
in the projection angle a. We compute the probabilities associ- 
ated to each solution assuming that the region of Msys values 
between uncertainties are equally probable for individual solu- 
tions: f-Bia ~ 83.5%, f Bib ~ 16.5%, Pug ~ 6 X 10-''%. 

Between the two possible bound solutions, a ~ 21 and 76 
degrees, the second one is rather unlikely since it means a dis- 
tance of ~ 3.7 Mpc between A115N and A115S, i.e. well 
larger than the virial radii, while the second solution means a 
distance of ~ 1.0 /i^q Mpc in agreement with a certain degree of 
interaction. When assuming a = 21 degrees the present collid- 
ing velocity is very large (~ 4400 km s"') and the cluster cores 
will cross after 0.08 Gyr. 

The characterization of the dynamics of A115 using these 
models is affected by several limitations. For instance, possible 
underestimates of the masses - e.g., if the subclusters extend out- 
side the virial radius - lead to binding probabilities larger than 
those computed above. The models do not take the mass distri- 
bution in the subclusters into account when the separation of the 
subclusters is comparable with their size (i.e. at small a) and do 
not consider the possible effect of small, low velocity groups. 
Moreover, the two-body model breaks down in a regime where 
Al 15N and Al 15S subclusters are already strongly interacting. 
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Fig. 15. System mass vs. projection angle for bound and un- 
bound solutions of the two-body model applied to A115N and 
A115S subclusters (solid and dotted curves, respectively, see 
text). The thin curve separates the bound and unbound regions 
according to the Newtonian criterion (above and below the 
curve, respectively). The horizontal lines give the observational 
values of the mass system and its uncertainties. 



Finally, the two-body model does not consider the possi- 
bility of an off-axis merger as suggested by the X-ray surface 
brightness distribution (Gutierrez & Krawczynski '2005) . For an 
analytical treatment which assumes that Al 15N and Al 15S are 
in a circular Keplerian orbit with the orbital plane perpendicular 
to the LOS, we refer to those authors. 

As for the low velocity groups, we consider a group corre- 
sponding to WGAPl likely centred around BCM-D and possi- 
bly interacting with the A115N-I-S complex. The values of the 
relevant observable quantities for the two-clumps system are 
then: the relative LOS velocity in the rest frame, Vr = 3376 
km s"' (as recovered from WGAPl and the mean velocity of 
CORE-A and CORE-B cluster rest frame); the projected linear 
distance between the two clumps, D = 0.89 /i^,} Mpc (as re- 
covered from BCM-D and the mean position between BCM-A 
and -B); the mass of the system obtained by adding the mass 
of a ~ 500 km s"' group to A115N-t-S, i.e. logMjys - 
15.4427;^ggg25- According to the Newtonian criterion for grav- 
itational binding we find a binding probability of 38%. The 
bimodal model indicates that the group is infalling onto the 
A115N+S complex with a merging axis intermediate between 
the LOS and the plane of sky (i.e., a ~40-60 degrees). 

5.2. Active galaxies 

It has been suggested that cluster-cluster collisions may trig- 
ger star formation in cluster galaxies (Bekki 119991 Moss & 
Whittle 2000; Girardi & Biviano 2002 and references therein). 
Caldwell & Rose d 19971 1 noticed that post-starburst galaxies 
are frequently found in clusters with evidence of past collision 



events. Here we analyse possible segregations between passive 
and active galaxy populations. 

Out of 85 cluster members we classify 48 galaxies finding 
34 "passive" galaxies ("k" in Table [TJ and 14 "active" galaxies 
(i.e., 9 "k-i-a"/"a-i-k" and 5 "e" galaxies, respectively), see also 
Fig.E] 

The velocity distribution of passive galaxies diff'ers from that 
of active galaxies at the 99.47% c.l. according to the KS-test, ac- 
tive galaxies having a larger mean velocity (see Table |2|. With 
present data this difference seems due to both Balmer absorp- 
tion lines ("k-i-a"/"a-i-k") and emission lines "e" galaxies. This 
suggests that the galaxy population of the high velocity subclus- 
ter (Al 15N) is very active with respect to other Al 15 galaxies. 
Indeed, we find that two out three classified galaxies in CORE- 
A are active galaxies (in particular BCM-A), while no active 
galaxy is found in CORE-B. That Al 15S has a rich, red galaxy 
population is also found by Rakos et al. (2000) suggesting that 
this system is evolutionarily "developed". Analysing the 2D dis- 
tributions we only find a marginal difference between passive 
and Balmer absorption lines galaxies (at the 92% c.l. accord- 
ing to the 2dKS-test) with three out of nine Balmer absorption 
lines galaxies closely located to BCM-A. Assuming that galaxy 
activity is connected with the merging phenomena, our results 
point out that A115N subcluster is more affected than A115S 
in agreement with our finding that A115N is less massive than 
A115S (see Sect.©. 

6. Summary & conclusions 

We present the results of the dynamical analysis of the rich, X- 
ray luminous, and hot cluster of galaxies A115, showing a bi- 
nary apparency (A115N and A115S) and containing a diffuse 
arc-shape radio emission, connected to A115N. This emission 
is considered a very anomalous relic since elongated relics are 
generally located at the cluster periphery, and with the largest di- 
mension roughly perpendicular to the cluster radial orientation. 

Our analysis is based on new redshift data for 115 galaxies, 
measured from spectra obtained at the TNG in a cluster region 
within a radius of 15'x 20' . We also use new photometric data 
obtained at the INT telescope in a field larger than 30'x30' . 

We select 85 cluster members around z ~ 0. 193 and compute 
a global LOS velocity dispersion of galaxies, cry - 1362^[Qg 
km s''. 

Our analysis confirms the presence of two structures of 
cluster-type well recognizable in the plane of the sky and shows 
that they differ by ~ 2000 km s ' in LOS velocity. The northern, 
high velocity subcluster (Al 15N) is likely centred on the second 
brightest cluster galaxy (BCM-A, coincident with radio source 
3C28) and the northern X-ray peak. The southern, low velocity 
subcluster (A115S) is likely centred on the first brightest clus- 
ter galaxy (BCM-B) and the southern X-ray peak. We estimate 
that Al 15S is slightly dynamically more important than Al 15N 
having a velocity dispersion of o-y = 900 - 1100 km s ' vs. 
(Tv = 750 - 850 km s"'. The virial mass estimates for the two 
subclusters are Mn(< Rvir.N = 2.2 - 2.5 /z^^Mpc) = 8.0-11.6 
xlO'^ MoandMs(< Rvk.s = 2.6-2.9 /i^^Mpc) = 13.8 
-18.9 xlO'4 Mo. 

Moreover, we find evidence for two small groups at low ve- 
locities. In fact, the galaxy distribution obtained from our INT 
data shows the presence of two concentrations (C and D) in 
addition to A115N and A115S (A and B). The C group sur- 
rounds the third brightest cluster member and coincides with 
the clump already found by Beers et al. (119831 see also our 
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Dressler-Schectman analysis). The D group surrounds another 
bright galaxy (BCM-D in this paper) and likely coincides with 
a faint X-ray emission shown by Chandra data (Gutierrez & 
Krawczynski 2005). 

Considering the complex structure of A115 we estimate a 
global cluster virial mass of 2.2-3.5 xlO^^ h^^ Mq. 

Our results agree with a pre-merging scenario where Al 15N 
and A115S are colliding with a LOS impact velocity Avrf ~ 
1600 km s ' . The most likely solution to the two-body prob- 
lem suggests that the merging axis lies at ~ 20 degrees from the 
plane of the sky and that the cores will cross after ~ 0. 1 Gyr. 

In our scenario where 1) A115S is more massive than 
A115N; 2) the dynamically important axis is the axis connect- 
ing the two merging subclusters; 3) this axis is perpendicular to 
the largest dimension of the relic, the anomaly of the All 5 ra- 
dio diffuse source is likely overcome supporting its relic nature 
(Govoni etal. |2001bl l. 
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Fig. 1. R-hand image of Al 15 taken with the WFC camera of the INT. Targets with successful velocity measurements are labeled 
as in Table [T] Circles and boxes indicate cluster members and non member galaxies, respectively. 



Table 1. This data will be availabe just after publication in A&A. 
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Fig. 2. /?-band image of the cluster Al 15 with, superimposed, the smoothed contour levels of the Chandra X-ray image ID#3233 
(photons in the energy range 0.5-5 keV; blue) and the contour levels of the NVSS (Condon et al. [1998 ) radio image (red). Letters 
A, B, C and D indicate the three brightest cluster members each one corresponding to Beers et al. ([1983) clumps and a fourth bright 
galaxy (see text). 
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Table 1. Continued. 
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